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Abstract 



We examine the effects of spin-orbit and spin-spin coupling on the inspiral 
of a coalescing binary system of spinning compact objects and on the gravita- 
tional radiation emitted therefrom. Using a formalism developed by Blanchet, 
Damour, and Iyer, we calculate the contributions due to the spins of the bod- 
ies to the symmetric trace-free radiative multipole moments which are used 
to calculate the waveform, energy loss, and angular momentum loss from the 
inspiralling binary. Using equations of motion which include terms due to 
spin-orbit and spin-spin coupling, we evolve the orbit of a coalescing binary 
and use the orbit to calculate the emitted gravitational waveform. We find 
the spins of the bodies affect the waveform in several ways: 1) The spin terms 
contribute to the orbital decay of the binary, and thus to the accumulated 
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phase of the gravitational waveform. 2) The spins cause the orbital plane to 

precess, which changes the orientation of the orbital plane with respect to an 

observer, thus causing the shape of the waveform to be modulated. 3) The 

spins contribute directly to the amplitude of the waveform. We discuss the 

size and importance of spin effects for the case of two coalescing neutron stars, 

and for the case of a neutron star orbiting a rapidly rotating WMq black hole. 
04.25.Nx, 04.30.-w, 04.80.Nn, 97.80.Fk 
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I. INTRODUCTION AND SUMMARY 



Coalescing binary systems of compact objects are the most promising source of gravita- 
tional waves which could be detected by laser interferometric gravitational wave detectors 
such as the currently funded U.S. LIGO and French/Italian VIRGO detectors These 
systems consist of neutron stars or black holes whose orbits decay because of the dissipative 
effect of gravitational radiation. The binary pulsar PSR 1913+16 is an example of such a 
system and has given us our first evidence that gravitational waves exist @. Laser inter- 
ferometric gravitational wave detectors will be able to observe the very late stages of the 
inspiral of coalescing binaries (typically the final several minutes) as the gravitational wave 
frequency sweeps through a detector's bandwidth from 10 Hz to 1000 Hz. 

In the previous papers in this series |3|-|6| we have studied the evolution of coalescing 
binary systems using a post-Newtonian approximation. The post-Newtonian approximation 
involves an expansion of corrections to Newtonian gravitational theory with an expansion 
parameter e ~ v 2 ~ m/r, which is assumed to be small (we use units in which G = c = 1), 
where m = mi + ni2 denotes the total mass of the system, and where r and v are the orbital 
separation and velocity. Since we are interested in the late stages of the inspiral, where 
the fields may not be so small, and the velocity not so slow, we use an expansion carried 
out to the highest practical order. We use equations of motion carried out to (post) 5 / 2 - 
Newtonian order, the order at which the dominant gravitational radiation-reaction damping 
forces occur. Schematically the equations of motion are given as 



where x = x x — x 2 denotes the separation vector between the bodies and r = |x|. We use 
a gravitational waveform carried out to (post) 3 / 2 -Newtonian order beyond the quadrupole 
formula. Schematically, 



d 2 x/dt 2 = -(mx/r 3 )[l + 0(e) + £>(e 3 / 2 ) + £>(e 2 ) 

+ (e 5 / 2 ) + ■■■}, 



(1.1) 





3 



where Q l i represents the usual quadrupole term (two time derivatives of the mass quadrupole 
moment tensor), D is the distance between the source and an observer, and TT denotes that 
the transverse traceless part of the tensor should be taken ((post) 2 -Newtonian contributions 
have recently been derived by Blanchet et al ]7|). Given a set of initial conditions, we evolve 
the orbit with the equations of motion and calculate the emitted gravitational radiation. In 
previous papers we ignored the effects on the motion and radiation due to the spins of the 
bodies. It is the purpose of this paper to take these effects into account. 

The contribution of the spins of the bodies to the equations of motion has been studied 
by numerous authors @-ll|. They include a contribution due to a spin-orbit interaction, 



and a contribution due to a spin-spin interaction. In Sec. |I| we write down the spin-orbit 
and spin-spin contributions to the equations of motion and show that, although they are 
formally post-Newtonian corrections to the equations of motion, for compact bodies they 
are effectively (post) 3 / 2 -Newtonian and (post) 2 -Newtonian corrections respectively, in part 
because they involve factors of order {R/r) where R is the size of the body which is on the 
order of m for a compact body. We add these contributions to our previous equations of 
motion, and obtain equations of motion valid for arbitrary masses and spins. Our equations 
of motion neglect tidal effects; for compact binary systems these effects are expected to be 
very small until the very late stages of inspiral p,|12|. We also ignore rotationally induced 
quadrupole effects which would enter at the same order as spin-spin effects; these also are 
expected to be small until very late stages, except possibly for very rapidly rotating Kerr 



black holes ||12|| . The major effect of the spins on the orbital evolution is that they cause the 
orbital plane to precess, thus changing its orientation in space. The spins themselves also 
precess; the precession equations are also written down in Sec. |T[ 

Recently, we calculated the spin contributions to the symmetric trace-free radiative mul- 
tipole moments which we used to calculate the spin contributions to the energy lost from a 
binary system due to the gravitational radiation emitted from the system |13|]. In Sec. |ITT] 
we present the details of that calculation, and calculate the spin contributions to the gravi- 
tational waveform, the angular momentum lost from the system, and the linear momentum 



ejected from the system. 

Since the emission of gravitational radiation tends to circularize the orbits of the binary 
system ||, we study the system under the assumption of circular orbits. This assumption 
should hold until the very late stages of inspiral when either the bodies start to merge, tidal 
effects become important, or the innermost stable circular orbit is reached in which case 
the system undergoes a transition from inspiral to plunge [pjl4|. Assuming circular orbits, 



we calculate the spin-orbit and spin-spin contributions to the energy loss and variation of 
the orbital frequency, and hence to corrections in the accumulated orbital phase. In Table 
| we compare these contributions to the analogous contributions due to the quadrupole 
term and other post-Newtonian terms for various binary systems. These contributions to 
the orbital phase are important since, in order to extract information from the observed 
waveforms, theoretical templates must match the observed waveform to within about one 
cycle (which corresponds to half an orbit) over the several hundred cycles that appear in the 
sensitive region of a detector's bandwidth (between roughly 40 Hz and 100 Hz). Note that 
the templates do not have to match the observed waveform within one cycle over the entire 
bandwidth from 10 Hz to 1000 Hz (as suggested by Cutler et al ]T5J) since the detector is 
more sensitive to the signal at some frequencies than in others [JL6|J . As seen in Fig. 2 of Ref. 



17| , 60% of the signal-to-noise from a binary inspiral will accumulate between 40 Hz and 
100 Hz due to the shape of the LIGO noise spectrum, and to the fact that there are more 
orbits at lower frequencies. In Table p] we list the orbital phase contributions for both the 
restricted bandwidth (40-100 Hz) and the entire bandwidth (10-1000 Hz) of a LIGO-type 
detector. We see that the spin-orbit correction to the accumulated orbital phase is important 
(unless the spins are very small), while the spin-spin contribution is negligible unless rapidly 
rotating black holes are present. Note that the contributions of the post-Newtonian terms 
are less important in the narrower bandwidth than in the entire bandwidth since the sensitive 
region of the bandwidth is at low frequencies, for which post- Newtonian corrections are less 
important in most systems. Conversely, the use of matched templates may allow estimations 
of spins via the spin-orbit terms |13| , [T5|] . 



Using the equations of motion and the equations of precession, we numerically evolve the 
orbit for a binary system of arbitrary masses and spins. Using the orbit we are then able 
to calculate the gravitational waveform emitted by the system. If the spins of the bodies 
are not aligned perpendicular to the orbital plane, the orbital plane will precess in space 
thus changing its orientation with respect to an observer. Since an observed gravitational 
waveform depends upon the orientation of the orbital plane with respect to the observer, 
this precession will cause the waveform to be modulated. Fig. |l| shows an example of this 
modulation. We see that these modulations depend significantly upon the observer's location 
with respect to the source, and upon the orientation of the detector relative to the incoming 
wave. The size of the modulations also depends upon the relative magnitudes of the orbital 
angular momentum L and the total spin S and upon their relative orientation. (See Fig. |2| 
for a description of the source coordinate system.) If the angle between L and S is small, 
the modulations of the waveform will be small. If |S| << |L|, then the modulations are 
small regardless of their relative orientation. For a circular orbit |L| = /i(rm) 1 '' 2 to leading 
order, while for each spin we define |Sa| = XA m \ where < xa < Xmax, where Xmax = 1 
for a black hole, and depends on the uncertain nuclear equation of state for neutron stars. 
For most neutron star models Xmax < 0.7 [IS|. In Fig. |3| we compare the relative sizes of 



L and S for the case of equal masses and for the case of a 10:1 mass ratio. We see that 
|S| << |L| almost always for the equal mass case, so that the modulations of the waveform 
will be small. For a neutron star orbiting a rapidly rotating massive black hole, however, the 
modulations can be substantial if L and S are sufficiently misaligned. The modulation of 
the waveform due to the spin-induced orbital precession has also been independently studied 
by Apostolatos et al |Tj|. 

There is also an explicit contribution to the waveform due to radiative multipole moments 
generated by the spins of the bodies. This contribution is relatively small until the late stages 
of inspiral. In Fig. ^ we compare the spin-orbit and spin-spin contributions to the waveform 
with the quadrupole contribution and higher-order post-Newtonian contributions for the last 
few orbits before coalescence. We see that the spin-orbit contribution is comparable to the 
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higher-order post-Newtonian contributions, while the spin-spin contribution is practically 
negligible. 

The rest of the paper presents the details. In Sec. |XT| we assemble the equations necessary 
to evolve the orbit. In Sec. [TTT] we derive the spin corrections to the radiative multipoles, and 
calculate the spin corrections to the gravitational radiation emitted, energy lost, angular 



momentum lost, and linear momentum ejected from the binary system. In Sec. |V| we 
examine the system in the limit of circular orbits, and in the limit of small precessions. 
In Sec. |V| we present our results for various numerically evolved orbits. Finally, in Sec. 



VT] we discuss these results, and their implications for extracting useful information from 
observations. In Appendix |A] we discuss the issue of "spin supplementary conditions" used 
to fix the center of mass of the spinning bodies to post-Newtonian order. In Appendix [FJ 
we list the post-Newtonian corrections to the circular orbit waveform. In Appendix we 
compare our results with calculations involving test masses orbiting spinning black holes. 

II. ORBITAL EVOLUTION EQUATIONS 

A. Equations of Motion 

Equations of motion for two bodies of arbitrary mass and spin have been developed by 
numerous authors (for reviews and references see [0-inj). By eliminating the center of mass 
of the system, we convert the two body equations of motion to a relative one-body equation 
of motion given by 

a = ayv + &pn + &so + &2PN + &ss + &rr, (2.1) 

where 

Tfl 

Sl N = — ^n, (2.2a) 

TTt f ' TTt 3 1 

SL PN = —={n {l + 3r ] )v 2 -2{2 + 7 ] )---r ] r 2 - 2(2 - V )fv \, (2.2b) 
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(2.2f) 



where ajv, &pn, and a 2 pjv are the Newtonian, (post) 1 -Newtonian, and (post) 2 -Newtonian 
contributions to the equations of motion, is the contribution to the equation of motion 
due to the radiation-reaction force, and &so and &ss are the spin-orbit and spin-spin contri- 
butions to the equations of motion which we have ignored previously, and where x = xi — x 2 , 
v = dx/dt, n = x/r, 77 = mim 2 /m, 77 = fx/m, 5m = 777,1 — ^2, S = Si + S 2 , and 
A = m(S 2 /m 2 — Si/mi), and an overdot denotes d/dt. 

It should be noted that the above expression for a^o is not unique; it depends on a 
"spin supplementary condition" (SSC) which is related to the definition of the center-of- 
mass world line x A for each body A. The above form of a 50 is for the covariant SSC given 
by Sa'uav = 0, where u A is the four-velocity of the center-of-mass world line of body A, 
and 



qflU 



X A 



X. 



(2.3) 



where r^ u denotes the stress-energy tensor of matter plus gravitational fields satisfying 
t m % = 0, and square brackets around indices denote antisymmetrization. Note that the 
spin vector S of each body is defined by S l A = ^€ijkS A k . We discuss the issue of SSCs in 



more detail in Appendix [A[ Here we simply wish to emphasize that since we have chosen a 
center-of-mass world line for each body through our choice of a SSC, we must ensure that 
all our calculations are consistent with this choice. 

Since the spin of each body is of order tjiRaVa where Ra is the size of body A and 
va is its rotational velocity, we see that the spin-orbit and spin-spin accelerations are of 
order (Ra/t)vva and (RA/r) 2 v\, respectively, compared to the Newtonian acceleration; 
these terms thus are formally of (post) """-Newtonian order. For compact objects, however, 
Ra is of order m, while va could be of order unity for sufficiently rapid rotation, so that 
the spin-orbit and spin-spin accelerations are effectively of (post) 3 ' 2 -Newtonian and (post) 2 - 
Newtonian order, respectively. If the bodies are slowly rotating, the spin contributions to 
the acceleration will be even smaller. 

It is interesting to note that while a^, &pn, &2PN, and slrr are all confined to the orbital 
plane, in general a^o and &ss are not. As a result, the orbital plane will precess in space 
(except for specific spin orientations) resulting in modulations of the observed waveform. 
We will discuss this effect in more detail in Sec. [TV| Iyer and Will J2t]] have derived 
post-Newtonian corrections to a^ at 0{e 7 ^ 2 ) and 0(e 4 ) where the latter are the spin-orbit 
corrections to radiation reaction. 



B. Spin Precession Equations 

In addition to the precession of the orbital plane, there are precessions of the spin vectors 



themselves. This effect has been studied by numerous authors |8HlO|; the relevent equations 
are 

S 1 = i{(L N xS 1 )(2+^)-S 2 xS 1 
r 6 I 2 nil 

+ 3(n-S a )nxSi), (2.4a) 



S 2 = -t((L n x S 2 )(2 + -^i) - Si x S 2 

+ 3(n- Si)n x S 2 1, (2.4b) 
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where Ln = /x(x x v) is the Newtonian orbital angular momentum, and where the first 
term in each expression is the precession due to spin-orbit coupling, while the second and 
third terms are due to spin-spin coupling. It is straightforward to show that the total spin 
S evolves as 



f7 a 35m " 

L N x -S H A 

1 2 2m, 



+ 3(n-Si)(n x S 2 ) 
+ 3(n.S a )(fixSiU. (2.5) 



It is useful to note that the precession of spins is a post-Newtonian effect, since L^/r 3 
(v/r)(/j,/r) Ri e(d/dt), and Si/r 3 ps rriiRv/r 3 pa e 3 ^ 2 (d/dt). 

Since the precession equations have the form Sa = x Sa, the magnitudes of the 
spins remain constant. The spin Sa instantaneously precesses about the vector with a 
precession frequency given by tu^ = \CIa\- It is an instantaneous precession since Qa is 
precessing itself in some complicated manner. Notice that if both bodies are spinning, the 
total spin S (with rare exceptions) does not have constant magnitude as the spins precess. 

C. Constants of the Motion 

Through (post) 2 -Newtonian order, the equations of motion can be derived from a gener- 
alized Lagrangian, that is a Lagrangian which is a function not just of the relative position 
and relative velocity, but also of the relative acceleration pi] . In our previous papers PJT3] 



we transformed the Lagrangian into relative coordinates and used it to compute the energy 
and total angular momentum of the system which are conserved to (post) 2 -Newtonian order, 
in the absence of radiation reaction. Combining our expressions for the non-spinning case 
and the spinning case, the energy is given by 

E = E N + E PN + E S o + E 2 pn + E S s, (2.6) 

where 

fir, ml 

E N = ii{-v 2 --j, (2.7a) 
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E PN = n{-(l- 3 V )v 4 + - (3 + i^t; 2 - + -17-r 2 + -(-) 2 j , (2.7b) 
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(2.7c) 
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+ |(4 + 69?7 + 12?7 2 ) ( 7 )^ 2 }, 

f? ss = ^{3(fi-S 1 )(n.S 2 )-(S 1 .S 2 )}, 
and the total angular momentum is given by 



mV 
r / 



J = L + S, 



(2.7d) 
(2.7e) 



(2.8) 



where 



L — ~L N + Lpjv + L50 + L 2 pjv, 



JN {|^ 2 (l- 377) + (3 + 77)^}, 
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V m 



Upn = L N { jj(l - 7 V + 1377 2 )^ 4 - ^(2 + 5??)™r 2 



+ -(7- 10n - 9n 2 )—v 2 
2 " r 

+ 1(14-41^ + 4^) 



(2.9a) 



(2.9b) 



;2.9c) 



(2.9d) 



Note that there is no spin-spin contribution to J. It is straightforward to show that to (post) 2 - 
Newtonian order, E = 3 = 0, where it is understood that whenever the relative acceleration 
is found in the time derivative, the equation of motion carried to the appropriate order is 
substituted. 
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D. Precession of the Orbital Angular Momentum 



Since the total angular momentum J is conserved (in the absence of gravitational radi- 
ation), it is clear that the orbital angular momentum L must precess as 

L = -S, (2.10) 

where S is given by Eq. (|2.5|) . 

If we restrict ourselves to the case of one spinning body then 

S = ^{i(l + 3— )(L N xS)|, (2.11) 

where m s is the mass of the spinning body. Since to lowest order J = L N + S, then 

§ = ^ (1 + 3 ^ )(JxS) }- (2 ' 12) 

Similarly, since L = Ln to lowest order, 

L = ^ (1 + 3 ^< JxL >}- (2 ' 13) 

These two equations imply that L and S precess about the fixed vector J at the same 
rate with a precession frequency given by 

^=|i(l + 3^). (2,4) 

Note that Ln is not necessarily parallel to L because of the Lgo terms, so that the orbital 
plane (determined by Ln) does not precess in the simple manner above. Instead the varying 
L S o terms cause it to wobble slightly on an orbital timescale as it precesses about J. This 
is illustrated in Fig. |. 

In the case of two spinning bodies, L, Si and S2 precess in a very complicated manner 
(with few exceptions), which can only be examined numerically. Fig. |^ shows an example 
of such a precession. Note that in this case the orbital plane tilts back and forth as it 
precesses about J. In Sec. [V|we will examine the precession of the spins and orbital angular 
momentum in the case of nearly circular orbits, and examine the effect of gravitational 
radiation on the simple precession for one spinning body, and on the more complicated case 
of two spinning bodies. 
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III. GRAVITATIONAL RADIATION EQUATIONS 



A. Symmetric Trace-Free Radiative Multipoles 

Our goal is to calculate the effects of the spins of the bodies on the gravitational radiation 
waveform emitted by the inspiralling binary and on the energy, angular momentum, and 
linear momentum radiated from the system. Thorne [23] showed that these quantities can 



be calculated using symmetric trace-free radiative multipole moments. For example the 
gravitational waveform to (post) 3 / 2 -Newtonian order is given by 



9 f ( 2 ) 1 (3) 1 ( 4 ) 

hv = — {pi +_ pi k N k + — P jkl N k N l 
D\ 3 12 



i (5) 

60 



_|_ e kl(i 



A (2) 1 (3) 

_ v) k N l + - ,V )km N l N r ' 
3 2 



2 (4) 

_j jj)kmn J^ljyrnjyn , _ _ 
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> , (3-1) 

TT 

where J u " are the mass multipole moments (see below), J 1 - 7 '" are the current multipole 
moments, D is the distance from the source to the observer, N l is a unit vector from the 
center of mass of the source to the observation point, the notation (n) over each multi- 
pole moment denotes the number of derivatives with respect to retarded time, e Jjfc is the 
completely antisymmetric Levi-Civita symbol, and parentheses around indices denote sym- 
metrization. To the accuracy we need, the mass quadrupole moment _P J needs to be calcu- 
lated to (post) 3 / 2 -Newtonian order beyond the lowest order, the mass octopole moment J yfe 
and current quadrupole moment J jfc need to be calculated to (post) ^Newtonian order, Pi kl 
and Ji km need to be calculated to (post) ^-Newtonian order, and pi klm and Ji kmn need to 
be known only to lowest order. 

Blanchet, Damour, and Iyer [p3] ,Z^[ (BDI) have developed a formalism for calculating 



these radiative multipole moments in terms of integrals over the source stress-energy. We 
use the BDI formalism to evaluate the spin-orbit and spin-spin corrections to the radiative 
multipole moments to the necessary order. 
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1. Mass multipole moments 



The mass multipole moments in harmonic coordinates (the coordinates in which our 
equations of motion are written) are given to post-Newtonian order by Eq. (3.34) in Blanchet 
and Damour as 



where L denotes a multi-index (i.e. x L = x %1 x t2 . . .x tl ), the superscript STF denotes that 
only the symmetric trace-free part is to be taken, and the source densities are given by 



In previous calculations involving nonspinning bodies we were able to evaluate the in- 
tegrals by assuming a point particle limit. Taking into account the spins of the bodies, 
however, precludes this. Instead we will assume the bodies to be well-separated, approxi- 
mately spherically symmetric (in harmonic coordinates), stationary, rigidly rotating compact 
objects whose structure is given by that of a perfect fluid. We will then neglect any effects 
due to the finite size of the bodies with the exception of each body's spin. The stress-energy 
tensor for a perfect fluid to the order we need is 




(3.2) 



a(x, t) = T 00 + T u , 



(3.3a) 



<r i (x,t)=T 



(3.3b) 



(3.4a) 



T 



(3.4b) 



T y = p*v l v 3 + pS lJ 



(3.4c) 
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where p* = p(l + \v 2 + 3U) is the so-called "conserved density" (it satisfies a continuity 
equation to post-Newtonian order) |[25|| , with p the local mass density, v the velocity, and U 
the Newtonian gravitational potential; II is the specific internal energy density, and p is the 
pressure. 

Substituting Eqs. ( |3.3| ) and ( |3^ ) into Eq. ( |3.2| ) for the case I = 2, the mass quadrupole 
moment is given by 



(iV) OT V(x) 



1 + ^ _ u + 3- 



P 



P W 



n 



j3 1 ^ 2 

o H — 

14 eft 2 



X X' 



x 2 p*(x)cf x 



20 d 



(x i x i x k ) STF p*(x)v*cPx 



(3.5) 



21 dtJA 

Following previous post-Newtonian calculations P,p6|, we choose the following provi- 
sional definition for the center of mass for each body: 



x A 



1 



x l p* (x) 



i + lv* A + n-±u A 



d 3 x, 



(3.6) 



where 



m A 



P W 



1 + -v\ + n - ^C/ A 



1 

— i 

2 



(3.7) 



where 



v l — v A , v\ = dx A /dt, and U A is the Newtonian potential produced by the 



A-th body itself. It turns out that this definition of the center-of-mass world line does not 
correspond to the center-of-mass world line of our equations of motion chosen through the 
use of a SSC, but rather to one given by a different SSC. There does exist a transformation 
between the two world lines, given by [E7I] 



x , 



x\ + 



1 



2m, 



(v A x s A y 



(3.8) 



This shift in the world line is of post-Newtonian order, so it can be neglected at lowest order. 
We choose to use our provisional definition of the center of mass to evaluate the integrals, 
and then use the transformation on the result so that it is consistent with our equations of 
motion. See Appendix |A| and Ref. p7| for more details. 
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Eq. ( p.5|) has been evaluated by several authors @,^| for the case of nonspinning bodies, 
i.e. v A = 0. By substituting x % = x l A + x\ and v % = v\ + v\ into Eq. ( |3.5| ), using the center- 
of-mass definition Eq. ( |3.6| ) and a virial theorem, and neglecting terms containing x A x 3 A 
which are 0(/3 2 ) relative to x l A x\ where (3 = R/r, Blanchet and Schafer |^6| have rewritten 
Eq. O as 



2* = £ \m A {x^ A ) STF 



B^A 



m B 
tab 



20 d 
21dt 



3 / p\x)v A {x i A x> A x k A ) STF d i x 



+ m A [x A x A x A 
1 



J ^fc \STF„,k 



+ 6v A lp*(x)(x A ^ A ) STF v\d 3 x+^ 



m A x A ( x A x a) 



(j/2 u p*{x)x A d 6 x + 2 J A P*(x){x\xtt 



l TJ 4- 3P 
2^ + 7W 



(fx • 



(3.9) 



The third integral in Eq. (|3.9| ) is just the mass quadrupole moment of body A, which 
is of 0(/3 2 ) so we will neglect it. The last integral in Eq. ( |3.9| ) will also vanish because 
of our assumptions of approximate spherical symmetry and rigid rotation. This just leaves 
integrals of the type: 



*/ \-a-b j3 / */ \ (-[a-b] . _(o_6)"l j3 

p ('x.)x A v A d x = J p (x) \ x A v A + x A v A j d x 
= \st + \jJ A P*{x)x A x h A d*x, 



where we have used Eq. fl2.3|) evaluated to lowest order. The second term on the right- 
hand side vanishes by our assumption of stationary spherical symmetry so we will neglect 
it. Changing the spin tensor to a spin vector we obtain 

1 



p*(x)x A v A d*x = -e^S\. 

A A 



(3.10) 



Substituting Eq. ( |3.10|) into Eq. ( |3.9|) and carefully counting the STF parts we obtain 



i J \STF 



J* = E \™a{x a x> a ) 

A I 

1 d 2 



1 + 



E 



B^A 



tab 



20 d 

21 ~dt 



m A(x A X A X A ) V A 



+ 



14 dt 2 



2 („A ^ STF 



m A x A (x A x 



+ 3 



x A (v A x S A ) j 



STF 4 d 
~ 3~dl 



^a( x a x S A ) j 



STF 



(3.11) 
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Using Eq. ( |3.8|) so that we have a consistent center-of-mass definition with our equations 
of motion, the mass quadrupole moment becomes 



i J \STF 



3 2 ^ m B 

2 B^A r ^ 



20 d r 

21 dt 



m A (x> A 4) STi? 4 



1^1 
+ 14 dt 2 



™a4(44) stf + 4 x s A 



stf 4 d 
" 3dt 



^a( x a x S A ) j 



STF 



(3.12) 



We rewrite the mass quadrupole moment in relative coordinates by using the transfor- 
mations 



Xl 



x 2 



m 2 : 1 5m / 2 m 



m 2 m 



v- j 

r 



mi 1 5m, 9 m N 
m 2 m r 



+ ^-(v x A) 

m 



+ -(vx A), 

m 



(3.13a) 



(3.13b) 



which can be obtained from a constant of the motion that can be taken as the center of 
mass 



2S 



29 . , 9 1 , , m 

l + ^(l-3^ 2 --(5-8^ 7 



The relative mass quadrupole moment to (post) 3 / 2 -Newtonian order is 

- ^(l-3ri)fxrr(x l v j ) STF 

+ 11(1 _ 3r])fir 2 (vV) STF + -7] [z* (v x ' 
2 X 3 

4 r , . . . .-l 5TF 



I TV 

1 J TajZ5 



(3.14a) 



where £ = S + —A, and where 7y aiZ is a reminder that tail effects need to be included at 
(post) 3 / 2 -Newtonian order. See Wiseman |J and Blanchet and Damour for more details 
on gravitational wave tails. 

Note that the spin-orbit correction is a (post) 3//2 -Newtonian correction for compact ob- 
jects. Since a spin-orbit contribution to a multipole requires a term involving va, and a 
spin-spin contribution vaVb, it is easy to see that Eq. ( |3~2] ) implies that there are no spin- 
orbit corrections to the higher mass multipole moments at lowest order, and no spin-spin 
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contributions at lowest order or at post-Newtonian order. Therefore the remaining mass 
multipole moments to the appropriate order are as given by Wiseman H 



jijk _ _ 



5m 



/'- 



m 



| (arVa; fc ) 



1 . , m 

_ - 5 - 1377 - 
6 r . 



STF r 1 . , 9 

1 + -(5-19?7)v 2 
6 

(1 - 2r ] )rr{x i x j v k ) STF 



+ (1 - 2r 1 )r 2 {x i v ] v k ) STF Y 



(3.14b) 



r jM = /i(l - 3r/) ( 



rr$ rfj rfM rfX 
%aj JU JU JU 



STF 



5m 



m v 



rv>3 tyJ 1 * rpl syTYl 
JU JU JU JU 



STF 



(3.14c) 
(3.14d) 



2. Current multipole moments 



The current multipole moments are given to lowest order as 

STF 



jiL = ^iab J a b x aL d 3V 



(3.15) 



At leading order, the current quadrupole moment J l] is a moment of angular momentum 
density; thus it will give an orbital angular momentum contribution as well as a spin con- 
tribution (effectively at 0{e 1 ^ 2 )). Although there is also a post-Newtonian correction to the 
current quadrupole moment, spin-orbit terms arising from this correction will be effectively 
a (post) 3//2 -Newtonian correction, which is higher than we need, since the contribution of 
J %3 to the waveform is already 0(e 1 ^ 2 ) at leading order. 

Substituting Eqs. ( |3.3j ) and (|3.4| ) into Eq. ( |3.15| ) we obtain 

-i STF 



J 13 



E 

L A 



P 



(x)v b x a x j d 3 x 



(3.16) 



for the current quadrupole moment. 



Substituting x l = x\ + x' l A and v l = v l A + v\ into Eq. (|3.16|) , using the center-of-mass 
definition Eq. (|3.6|) , and neglecting terms of 0(j3 2 ), we obtain 



I A 



m A x A x A v A + x A I p*(x)x A v A d x 

J A 



+ x A I p*(-x)x A v b A d 3 x 



STF 



(3.17) 



Using Eq. (|3.10|) the current quadrupole moment is then given by 

.i STF 3 



x A (x A x v A ; 



+ ^ ( x aS J a) 



STF 



(3.18) 



Note that the spin-orbit correction is effectively (post) 1/,2 -Newtonian order. Repeating the 
calculation for the current octopole moment we obtain 



x A x A (x A x v A ) fc 



STF 



2 (x l A x A S\ 



STF} 



(3.19) 



Transforming Eqs. (3.18) and ( |3.19| ) into relative coordinates and adding the post- 
Newtonian correction to the current quadrupole moment derived by Wiseman Q, we obtain 



T 3 



5m 



-H- 



X X X V 



m I 

I Tfi 

+ —(27 + 30r7) — 



STF 



l + l(13-68^ 2 



x 



14 

?/(x x v) 



r 

STF 



+ ^( 1 -2r?)rf 



(3.20a) 



J^ k = n(l-3 v ) 



a;V(x x v) fc l STF + 27/(xV«£ fc ) STF 



(3.20b) 



The final moment we need is given by 



Jiiu= 5m 
m 



X l X j X k (TL X vV 



STF 



(3.20c) 



Note that there are no spin-spin contributions to the current multipole moments at this 
order. 
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B. Gravitational Waveform 



Taking time derivatives of the radiative multipole moments ( 3.14Q and (|3.20|) , substitut- 
ing the equations of motion ( |2.1|) where appropriate, and substituting the results into Eq. 
the gravitational waveform is given by 



h 13 



2/i 



Q ij + p 0.5 Q ij + pQ ij + p( yj (j + pl.5 Qlj 



+ P™Q%o + P'-'Qtu + p 2 Q ij 



ss 



TT 



(3.21) 



where 



v v J n n J 



Sm (m 



p0.5 Q i j = ^ 3^ 2n (i vj ) 



7U 



rn l nP 



(N • n) + 



—n l n j — 2v l v j 
r 



PQ» = -(l-ty) 4 
3 I r 



m 



3rn l n j — 8n^v J ^ 



(N-n)(N- v) + 2 



(3.22a) 

(N-v)J>, (3.22b) 
(N • v) 2 



+- 



+ 



m 



(3v 2 - 15r 2 + 7— )nW + 30rn {i v j) - 14uV 
r 



3t>V — —n l n j 
r 

(N • n) 2 l + -— f (5 + 3r])n {i v j) 
) 3 r 



m 



(i-3 V y- -(2-377) 

3 r 



vv J + 



m 



(l-3r ] )r 2 - 1 -(10 + 3 V y + f^ 



n % n 3 , 



PQlo = -(A x N)fcn- 



«rj.J') 



(3.22c) 
(3.22d) 



pl . 5 Q« = ^ (1 _ 2 ,){i- 
m 4 r 



(45r 2 - 9v 2 - 28— )nV + 58vV - 108rnM 
r 



(N-n) 2 (N- v) 



1 

+ 2 



— wW - 4t/t> J 
r 



(N-v) d + - 

r 



+ 44— )n ( V } - — fuV 
r ' 2 

5m 1 m ~ A J : ,■, 

H (N • n)< nWr 

m 12 r v y I 



5 777- I 

-(3v 2 - 7r 2 + 6— )rn i n j - -(21w 2 - 105r 2 
4 V r' 6 V 

3 m 



(N-n) 3 + lOnM - 3rnV' (N-n)(N-v 



2 r 



m 



r 2 (15 - 90?]) - w 2 (63 - 54r/) + —(242 - 24r/) 



-rv l v J 



(186 + 24r?) + 2n ( V ) 



m 



r 2 (63 + 5477) (128 - 36r/) + t; 2 (33 - I877) 



— (N • vK -vV 

m 2 



rn 



(3 - 8r?) - 2w 2 (l - 577) 



n {i v j) — f(7 + 477) 



r 



3 I 777- I 

(1 _ 2 77)r 2 + -(26 - 3r/)- - -(7 - 2tt> 2 

4 3 r 4 



(3.22e) 
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(n x v)-(12S + 6— A) 

m 



vx(9S + 5— A) 

m 



j) 



+ 



3f(N • n) - 2(N • v) 



+ frv 



, „ 5m , 
fix 12S + 6 — A 

m 



;S + ^A)xN 

m 

i) 



nx(2S + 2— A) 

m 



i) 



2(N • nl 



!S + ^A)xN 

m 



J) 



m 



/X •/ 



oo ( 4 ) 



V2s/ 12 



du' 



(3.22f) 



(3.22g) 



P 2 Qss 



{nV [(Si ■ S 2 ) - 5(n • Si)(n ■ S 2 )] + 2r^' ) (n ■ S 2 ) + 2n^' ) (n • Si)} , 

(3.22h) 



where Q y is just the standard quadrupole term, P°- 5 QV and PQ % i were derived by Wagoner 
and Will |(J, P L5 Q^ was derived by Wiseman f§, and PQ% and P l ' 5 Q% are the explicit 
spin-orbit corrections to the waveform. P^QtoU * s ^ ne leading order contribution to the 

■ ' STF 

waveform due to the tail, where = . Notice that it depends upon the past 

history of the binary's inspiral (u is a retarded time), and that s is an arbitrary matching 
parameter. See Wiseman || and Blanchet and Damour |2{| for more details about gravi- 
tational wave tails. Note that we have included the leading order spin-spin contribution to 
the waveform even though it is effectively a (post) 2 -Newtonian term. It is due entirely to 
substituting a S s into the time derivatives of the mass quadrupole moment. Note that we 
have simplified Eq. Q3.21 ) by using relations such as 



i j) b k 



STF 



N 



TT 



b x N) (l a 



TT 



Fig. |] shows the different contributions to the waveform for the final few orbits of an 
inspiralling binary system. We see that the spin-orbit contribution is comparable to the 
other post-Newtonian contributions, but the spin-spin contribution is almost negligible. 



C. Energy Loss 



The radiative energy loss in terms of STF radiative multipoles is given by Thorne |22 



as 
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dE _l/t 3 )® 
dt~ 5 lj!J + 



5 W (4) 16 (3) (3) 1 



(3.23) 



for the accuracy we require. Taking time derivatives of the radiative multipole moments 
G3.14T ) and ( |3.20| ), and substituting the equations of motion ( [2.1|) where appropriate, the 
energy loss is given by 



where 



dE 

— = E N + E PN + E S o + E TaU + E S s, 
dt 



E 



N 



m 2 fj? 



15 r 4 



{l2v 2 - llr 2 } 



2 m 2 /j, 2 
105 r 4 



E 



PN 



\ (785 - 852T?)?; 4 - 160(17 - r])—v 2 
I r 



m 



2-2 



+ 8(367 - 1577) — r 2 - 2(1487 - I392rj)v 2 t 



+ 3(687 - 620r/)r 4 + 16(1 - 477) ( 



(3.24) 



(3.25a) 



(3.25b) 



E 



so 



15 r 6 



S(78r 2 - 80v 2 - 8—) 



+ ^ A (51r 2 -43t; 2 + 4- 

m r 



(3.25c) 



2 (3) d 



(3.25d) 



,95' 



15 r b I v ' 

+ 3(Si • S 2 ) (47v 2 - 55r 2 ) + 71 (v • Si)(v • S 2 ) 

- 171r [(v • Si)(n • S 2 ) + (n • Sx)(v • S 2 )]}, 



(3.25e) 



where Ep^ was found by Wagoner and Will , and Eso and Ess were reported in our 
previous paper [[□[] ■ -Sraii depends upon the past history of the system and can only be 
evaluated explicitly for simple cases. Notice that the spin-spin contribution to the energy 
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loss, which is effectively a (post) 2 -Newtonian correction, comes from using post-Newtonian 
equations of motion in the derivatives of the mass quadrupole, and also from the contraction 
of the current quadrupoles. We have ignored (spin) 2 -terms which are the same order as the 
spin-spin terms. In Fig. [7| we compare the spin contributions to the energy loss with the 
other contributions for an inspiralling binary system. Again we see that the spin-orbit 
contribution can be significant, while the spin-spin contribution is almost negligible. 



D. Angular Momentum Loss 



The radiative angular momentum loss in terms of STF radiative multipoles is given by 
Thorne 1231 as 



d£_ 
~~dT 



2 f (2) (3) 5 (3) (4) 16 (2) (3) ~ 

— ~ e ijk j Ijilkl ~^~Y26 ^ lm ^ klm ~*~~g" ^fl^kl 



(3.26) 



for the accuracy we require. Taking time derivatives of the radiative multipole moments 
Q3.14j ) and (|3.20| ), and substituting the equations of motion ( j2.1|) where appropriate, the 
angular momentum loss is given by 



dJ • • • • • 

— — Jn + JPN + JSO + JTail + JSS, 

at 



(3.27) 



where 



1 N 



5 r 5 



2v z - 3f 2 + 2 



m 



(3.28a) 



\ PN = — — ^L N ((307 - 5A8r])v 4 - 6(74 - 277^)t» 2 r 2 + 2(372 + 197t?) — r 2 
105 r b I ' r 

+ 15(19 - 72ry)r 4 - 4(58 + 95r/)— v 2 - 2(745 - 2rj)(— ) 2 L (3.28b) 

r r J 



• 4/x 2 (2m,. 2 2N „... . .... „ 

J.so = — ^ < ir - v z ) A - r— nx 



5r 3 r 



5m 
m 



. m 
r 



., _,, 5 <5m. . , 

4(v x S) H (v x A) 

3 m 



m 
H nx 

r 



-rvx 



„. , 9 41 9 4m N 5m, AW 9 , 2m ( 
n x S 15r 2 - — v 2 + -—) + — n x A)(9r - 8v 2 - -—) 
6 6 r m 3 r 

x S)(18— + Uv 2 - 55r 2 ) + 5 — (n x A)( + Av 2 - 5r 2 

r m 3 r 
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-vx 



, . 9 71 9 50 m. 5m. AW 9 35 , m, 

(v x S)(36r 2 v 2 ) + — (v x A)(18f 2 v 2 - 9— 

3 3 r m 3 r ' 



Ln T 



(65r 2 - 37, 2 - ^™)S + (35r 2 - I9v 2 - ZI™)^ A 
3 r 3 r m 



(3.28c) 



(3) 



(2) 



(3.28d) 



Js5 = 



o |'(2) (3) (2) (3) 



w lfi (2) (3) 

a SS )+ /g (ass) J# (a N ) + - J&, (a N ) J« G (a N ) } , (3.28e) 



where Jpjy was calculated by Junker and Schafer pi| . Note that J Tan depends upon the past 
history of the system. To avoid lengthy expressions, we have left J$s m terms of derivatives 
of the multipole moments, where we have specified which contribution to the equations of 
motion should be substituted for the accelerations which appear in the time derivatives. 
Note that while Jjv and Jpat are in the direction of Ln, in general Jso and Jss are not. 



E. Linear Momentum Loss 



The radiative linear momentum loss in terms of STF radiative multipoles is given by 
Thorne [22] as 

dP i 



dt 



2 W (3) 16 (3) (3) | 



(3.29) 



for the accuracy we require. Taking time derivatives of the radiative multipole moments 
( |3.14| ) and ( |3.20| ), and substituting the equations of motion Q2.ip where appropriate, the 
linear momentum loss is given by 

dP 



dt 



Pn + Pso, 



(3.30) 



where 



N 



8 5m 2 f m\ 4 . . „ 
-i] — rn 



105 m \ r 



hhv 2 - 45r 2 + 12— 
r . 



+ v 



38r 2 - 50v 2 - 8- 



m 



(3.31a) 
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Pso = -^^{4r(v x A) - 2v 2 (h x A) 
15 r° L 

-(fix v) [3r(n- A) + 2(v • A)]}, (3.31b) 



where Pat was studied by Fitchett |32|. Note that Pso, which is effectively a (post) 1 / 2 - 
Newtonian correction, can be directed out of the orbital plane. There is no spin-spin con- 
tribution to the linear momentum loss at this order. Wiseman [Q has derived the post- 
Newtonian corrections to P. In Fig. [8] we plot the momentum ejected for a typical inspiral. 
Notice that even in the presence of spinning bodies, the momentum ejection is periodic, and 
therefore there is no large buildup of momentum ejected in a specific direction. 



IV. CIRCULAR ORBITS 
A. Circular Orbit Limit 

Gravitational radiation tends to circularize the orbit of an inspiralling binary. Therefore 
we would like to examine the last several minutes of the inspiral with the assumption that 
the orbit is quasi-circular, that is, the orbit is circular on an orbital timescale, but inspirals 
on a radiation-reaction timescale. This is a reasonable assumption, since Lincoln and Will 
|| have shown that virtually all captured binaries will have sufficient time to circularize 
their orbits before plunging to coalescence. 

The equations of motion can be rewritten using the identities 

fi • a = r — rco 2 , (4.1a) 



A • a = ru + 2rcu, (4.1b) 

L N -a= -ruA^V (4.1c) 

where A = Ln x n, Ln = Ln/|Ln|, and the orbital angular velocity u is defined by v = 
rh+ruX. A circular orbit on a fixed plane is given by the solution r — f — Co — (g?Ln / dt) = 0. 
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This solution exists if ru 2 = — n ■ a, A ■ a = 0, and Ln ■ a = 0, where we have substituted 
r = into the righthand sides of Eqs. (|4.1|) . Examining the equations of motion ( |2.1| ) (while 
ignoring the radiation- reaction terms), we see that circular orbit solutions exist only if the 
spins are aligned perpendicular to the orbital plane. If we instead define a circular orbit as 
one that has a constant orbital separation, but allow the orbital plane to precess, there exist 
circular orbit solutions for the case of one spinning body, but not for the case of both bodies 
spinning with general orientations. 

If we first average the spin-orbit and spin-spin terms in the acceleration over an orbit, 
then we can obtain orbits of constant separation for arbitrary spins and orientations. In 
order to average over an orbit we need to assume that the spins and orbital plane remain 
constant over an orbital period, in other words, that uj p /u) is small. For circular orbits 

/ m \3/2 

muJ = i—J , (4.2) 



to leading order. Using Eq. fl2.14|) for the case of one spinning body we see that 

OJr. 



Since |J| < |L| + |S|, \L\/m 2 
see that 



/ m 


\3/2 




h3^" 






1 - 




\ r 








i/ 2) 


and 


|S| 


jm 2 = 



(4.3) 



uj l v r 



m — 
77 + 3- 



m 



' m \3/2 m ^ 

+ X' 

r / m 



m s 
3 + — 

m 



(4.4) 

which is small (~ 0(m/r)) until the very late stages of the inspiral, where the whole circular 
orbit approximation breaks down anyway. We would expect a similar argument to hold for 
the case where both bodies are spinning, since the spins' instantaneous precessions have a 
form similar to that of the precession in the single spin case (see Sec. [11 B| ). An examination 
of numerical evolutions of the precession equations confirms this. 

In examining circular orbits, we will assume an orbit where f = and ru 2 = — (n • a) 
where the brackets denote an average over an orbit. We therefore obtain the following 
expressions for a circular orbit. The orbital velocity is given by v = ru where 
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2 2 

r uj 



(7)H-'>(7)-& 



- 2^X1X2 



x ,(L N -Si)(2^ + 377) 



m \ 3/2 

_ + 

r / 



m\ 2 
r J 



(§1 • s 2 ) — 3(L N • Si)(L N • s 2 
where Sa = Xa^a^a- The energy and angular momentum are given by 



41 9N 
(6+— rj + rj 2 ) 



(4.5) 



^ m ^3/2 



1 1 r ~ 

-(7 - 49?7 - r? 2 ) - -77x1X2 [(§1 • § 2 ) - 3(L N • §i)(L N • s 2 ) 



(7)' 



J = /i(mr) 1/2 L N 1 + 2 



ro\ 1 



41 i=l,2 L 



'"7 



Xi (L N .Si)(8^ + 777) 



r 77x1X2 



(§1 • s 2 ) - 3(L N • §i)(L N • s 2 ) 



+ S - >(mr)^ £ 

4 i=l,2 L 



X*si(4— I + 77) 

77T 



rn 
r 

/ m ^3/2 

V r 



77T 
2' 



m \ 3/2 
_ + 

r / 



(4.6) 
(5 - 97?) 



For a circular orbit, the waveform is given by 



2ji (m 



hl3 = T\-)\ Q " + P Q < [7 ) + PQ " [7 ) + P Q < [7 ) 



m \ 3/2 I 



TT 



where 



(4.7) 



(4.8) 



Q 



v — 9 
c — z 



po.tQij == . fi ) n (i A i) + (N • A) [nV - 2A l A J ] } , 



(4.9a) 
(4.9b) 



PQ« = -(1 - 377) (N-n) 



5nW - 7A l A J 



16(N-n)(N-A> (i A j) 



P^Qv = ^ j(i _ 2?7) i(N • A) 3 (raW - 4A*A j ) + ^(N • n) 2 (N • A) (58A l A J - 37raW') 



3A*A J - raV 



1 



+ -(19 - 3ti) (nV - A l A J ) + — jn^A x N) J '\ (4.9c) 



to' 



1 



— (N • n)V*A j) + 15(N • n)(N • A)V*A j) 



--(19 - 477)A^ 



(N-A) 

- -(149 - 677) (N • n)n (i A j) l - Al^A j 
6 ) m 2 \ 



— (101 - 12?7)nV 



L N -(5S + 3— A) 

to 



6n l n J 



L N .(2S + — A) 

to 



+ 2(N • A) 



(S + — A xN 

to 



+ 2A ( * 



- /r, $ m a X 

fix S + — A) 

m 



3) 



Ax(9S + 5— A) 
m 



V') 



n 3) + 2(N • n) 



(S + ^A)xN 

to 



X j) \, (4.9d) 
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where X 1 = v l /\v\ for a circular orbit. The energy loss and angular momentum loss for a 
circular orbit are given by 



32 2 / m 
dt 5 \ r 



1 



^(2927 + 420,)^) 



— y. 

12 ^ 

lz i=l,2 



9 



^(^•§0(73^ + 75^) 



Air 



m \ 3 / 2 
r 



223(§! ■ s 2 ) - 649(L N ■ §i)(L N • s 2 ) 



(4.10) 



dJ_ 
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+ 527/) 







) 4 (mr) 1/2 < 


Ln 



1 - ^(2423 + 588,,) (?) - Q £ 



r / 



1 



+ 24 S 

^ i=l,2 . 



X i s l (37-|+427 ? ) 



77T 



m s 3/2 

r / 



Xi (L N -Si)(53-^ 



(4.11) 



where the 4ir terms are due to the gravitational wave tail. 

The rate of inspiral is given by r = (dE/dt)/(dE/dr). Taking Eq. fl£lQp and dividing it 
by the derivative of Eq. ( |4.6|) , we obtain 

i2 



dr 64/m\ 3 [ 1 , , fm 

V[—) 1 1 - ^(1751 + 58877) 



dt 



5 

47T 



336 



- V 

12 ^ 



^(^•§0(19^ + 1577) 



3/2 



:TOiX2 



59(§i • s 2 ) - 173(L N • §i)(L N ■ s 2 ) 



2 N 



(4.12) 



r) 46 

The above expressions can be inverted to express everything in terms of the orbital 
frequency. For a given orbital frequency, the separation is given by 



(r/m) = (muo)- 2/ A 1 - -(3 - ?))(mu) 2/3 - - £ 
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Using Eqs. ( |4.12j ) and ( |4.13| ), we find that the evolution of the orbital frequency is given by 
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The evolution of the orbital frequency can be used to calculate the accumulated orbital 
phase of the binary. The orbital phase as observed by a phase sensitive detector such as 
LIGO/ VIRGO is given by 
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(4.15) 



where ti is the time at which the signal enters the sensitive bandwidth (corresponding to a 
lower frequency u>i set by seismic noise) and tf is the time at which the signal leaves the 
sensitive bandwidth (corresponding to an upper frequency ujf set by photon shot noise), the 
time at which the orbit begins to plunge (corresponding to a frequency Uf of the innermost 
stable circular orbit), or the time when the two bodies begin to coalesce. The result is 
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The term with 107r is the tail contribution to the orbital phase. Table | shows the contribu- 
tion of each term to the orbital phase for several cases. 

B. Radiation- Reaction Effects on the Precessions 

When averaged over a circular orbit, the orbital angular momentum and spins precess 



as 
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Since (Lso) = over an orbit, the precession of Ln is also given by Eq. ( |4.17a|) with L 
replaced by Ln. Thus we see that for a circular orbit the magnitudes of the vectors L, Ln, 
and Si are conserved on average during precession in the absence of radiation reaction. 

For a binary system inspiralling due to gravitational radiation, there is a loss of total 
angular momentum given by J. We assume that the individual spinning bodies are suffi- 
ciently axisymmetric that they will emit negligible gravitational radiation on their own, so 
that their spins are unaffected by radiation damping to the order which we are considering. 
(See Ref. (19| for a more rigorous argument.) This means that the angular momentum loss 
is entirely from the orbital angular momentum, and that the total change in orbital angular 
momentum is the sum of the precession and the radiation damping. 

Let us examine this in more detail for the case of one spinning body, where we only 
consider the leading order damping effects. Then we have (see Eqs. (|2.12|) - (|2.14|) ) 

S = u p J x S, (4.18a) 
L N = UpJ x L N - e m L N , (4.18b) 

j = -e RR L N , (4.18c) 

where J = J/|J|, and where err = (32/5)(yu/m 2 )(m/r) 4 is the rate of angular momentum 
loss due to gravitational radiation. Defining L = |Ln|, S = |S|, and k = cos" 1 (Ln • S^, 

J = | j | = |l 2 + S 2 + 2LS cos k} 1/2 , (4.19a) 

i = cos" 1 (L N • J) = cos" 1 {[L + Scosk] /J} , (4.19b) 
we see that Eqs. ( |4.18| ) imply that 

5 = 0, (4.20a) 



L — — £rrL, 
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(4.20b) 



J = — €rrL cos i, (4.20c) 

k = 0, (4.20d) 

di €rrL . . , . 

— = — — smi. 4.20e 
dt J y 1 

As the binary inspirals, Ln and J shrink and the angle between them grows. This is shown 
in Fig. |9]. The angle k between Ln and S remains fixed as they precess about J, so that 
S tips toward J and at late times when J « S, i approaches k. Working to first order in 
(L6rr)/(Juj p ), we can show that, if J Q is the initial direction of J, 



(Jo x L T 



£rr 



1 — (t — t ) COS I 



Jo, (4.21) 



valid for (t — t a ) << J/(LeRji), which implies that J spirals about J Q as it shrinks. As 
long as the ratio A = {Lerr) / (Ju p ) « 1, then J will remain relatively fixed in space, 
precessing on a tight spiral around its earlier direction. If the ratio A is not small, however, 
Ln and S will still precess about J, but J will start to tumble. An example of this would 
be if L and S were nearly equal in magnitude, but pointing in opposite directions so that 
J=|L + S|<<L. This is the case of "transitional precession" described in Ref. |L9 . 



If both bodies are spinning, the precessions of L and S are more complicated (as in Sec. 
Il|), but the effects of radiation damping are qualitatively the same. As the binary inspirals, 
Ln and J shrink and the maximum angle between them grows. J remains relatively fixed in 
direction as long as the precession timescale is shorter than the inspiral timescale, and as long 
as the ratio A is small, where we replace the precession frequency uj p with the instantaneous 
precession frequencies u p A \ Fig. [K| shows an example of this. Numerical evolutions of the 
precession equations in the presence of radiation-reaction agree with our description above. 



C. Wave Polarization States 

The gravitational radiation emitted by the binary can be written in terms of its two 
polarization states h + and h x . The polarization states can be expressed as linear combina- 
tions of the components of h 1 ^ in some suitable coordinate system. Normally one chooses 
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a coordinate system such that the orbital plane lies in the x-y plane, and the source and 
detector are located in the x-z plane as this simplifies the equations involved (see Ref. [|33|j). 
If the bodies are spinning, however, this coordinate system is not fixed in space. This does 
not prevent one from using such a coordinate system (see Apostolatos et al JT9|), but one 
must remember that it is precessing. Instead, we choose our coordinate system such that 
our z-axis lies along some initial direction of J. We define O = cos -1 • J Q ^ and choose a 
coordinate system in which z = J Q and N = cos Gz + sin Gx lies in the x-z plane (see Fig. 
HD- (Note that if G = 0, x can be chosen arbitrarily.) Following the method of Finn and 
Chernoff ]33| we define the radiation coordinate system such that 



ef = N, (4.22a) 



{1-(J„.NP} 



e ^ x N ^< J °- N > N - J °, (4,2c) 



{l-(Jo-N) 2 } 



Then from Eqs. (3.2) and (3.3) of Ref. [Q, the radiation can be written in terms of its 
polarization states 

h+ = \ {cos 2 Qh xx - h yy + sin 2 Bh zz - sin 20h xz ] , (4.23a) 

h x = cos Qh xy - sin Qh yz . (4.23b) 

The response of a detector will be the linear combination 

h = F + h + + F x h x , (4.24) 

where F + and F x are the antenna patterns which depend upon the orientation of the detector 
with respect to the binary. 

F + = |(1 + cos 2 9) cos 2(p cos 2^- cos 6 sin 2<f) sin 2ip, (4.25a) 
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F x = |(1 + cos 2 6) cos 20 sin 2^ + cos sin 20 cos 2^, (4.25b) 

where (#, 0) is the location of the binary with respect to the detector (whose arms lie along 
the x and y axes in the detector's coordinate system, with the z axis in the vertical direction), 
and ijj is the polarization angle between the gravitational waves' polarization axes and the 
direction of constant azimuth. (See Fig. [1]].) For our definition of the polarization axes, the 
polarization angle is given by 

V> = tan- ( 5 *<?° X }\ 5 ) . (4.26) 

Vjo-z-(j -n)(z-n); 

If N = ±z then 

y • Jo 



ip = tan 



-i 



-(N-z) 



(4.27) 



X ■ J 

If we define $ as the orbital phase with respect to the line of ascending nodes (the point at 
which the orbit crosses the x-y plane from below), then for a circular orbit the polarization 
states will be given by 

2/i f m\ o05n ( m\ 1 / 2 i n ^ ( m 



+ P ^ Q+ (^y"Y (4.28) 
where the quadrupole term is given by 

Q+ = -2[C+cos2$ + S + sin2$], (4.29) 
and similarly for h x where + is replaced by x, and where 



r~i 1 2 r\ I ■ 2 2- 2 \ , 1 / 2 • • 2 

G + = 2 cos B I sm a — cos tcos a I + ~ (cos zsm a 
— cos 2 a) — | sin 2 9 sin 2 z — | sin 26 sin 2i cos a, 

(4.30a) 

«S+ = | f 1 + cos 2 0J cos i sin 2a + | sin 26 sin % sin a, (4.30b) 
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C x = — I cos 6 sin 2a (l + cos 2 1) — ~ sin 6 sin 2z sin a, (4.30c) 

S x = — cos cos i cos 2a — sin 9 sin % cos a, (4.30d) 

where (i, a) are the spherical coordinates describing the direction of Ln (see Fig. |^). We 
give the post-Newtonian corrections to h + and h x in Appendix [FJ. The evolution of i and 
alpha is given by the precession equations ( [4.1 7|) . The evolution of $ is given by 



$ = co>-dcosi. (4-31) 

Note that our expressions for the polarization states (specifically the quadrupole terms) 
are much more complicated than the expressions in Apostolatos et al fl9|| because we have 



defined them with respect to a fixed coordinate system as opposed to a rotating one. In 
our description a varies as ~ u p t as the orbit precesses in the case of simple precession. On 
the other hand, in our description the polarization angle ip (|4.26|) (and thus the antenna 
patterns F + and F x ) is constant during the binary's inspiral, while in the description of 
Ref. [|19| it is not. The two descriptions are of course equivalent; we have simply made the 



complexity of the waveforms more explicit. 

From the above equations, it can be seen that the signal in the detector due to the 
quadrupole waveform can be written as 

h = C Q cos2<l> + SQsin2<l>, (4.32) 

where 

4/1 (m s 



C Q 1 n \ 

D \ r 



\ [C + F + + C X F X ], (4.33a) 



D \r 

The signal can be rewritten as 



Sq = (-) [S + F + + S X F X ] . (4.33b) 



h = A Q cos [2$ - S Q ] , (4.34) 
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where 



A Q = {C 2 Q + S 2 Q } 1, \ (4.35a) 
6 Q = tan" 1 (S Q /C Q ) . (4.35b) 



D. Small Inclination Angles 



Let us examine the waveform in the limit of a small inclination angle i. This limit will 
be valid if the total spin S is nearly aligned with the orbital angular momentum L or if 



|S| << |L|. In Fig. [12] we show the region for which the precession angle i is less than 
0.2 for an equal mass system. If we expand Eq. (|4.29 ) through 0(i 2 ), the overall signal 
amplitude due to the quadrupole term ( |4.35a|) will be given by 
2\i fm s 
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+ 4F 2 (9, 4>, ip) cos 2 6 + i sin 26 cos a 
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+ F + (9, 0, ip)F x (9, <f>, i>) -2% sin d 6 sin a + 3r sin 2 6 cos 6 sin 2a 



(4.36) 



Notice that to lowest order in i, the amplitude is constant, independent of the precession 
angle a, and that the modulations are of 0(i) and have the same frequency as the precession 
frequency, da/dt. For some detector orientations (see below), the 0{i) terms are suppressed 
and the modulations will be of 0(i 2 ) and go as twice the precession frequency. 

We can use Eq. ( [4.36| ) to explain the features of Fig. |I[ Note that this case has been 
presented in Cutler et al fl5]| , and in Apostolatos et al [|HJ for a single detector orientation 
which corresponds to 7 = i. In Fig. [I], 6 = 7r/2, F + = cos 27, F x = sin 27, and initially 
i w 0.1, so that 
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Q 



I cos 2 27 



2\i m 

— 2% sin 47 sin a + Ai 2 sin 2 27 j 
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1 — i 2 + 3i 2 cos 2a 
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(4.37) 



For 7 = 0, the modulations of the amplitude will have a frequency of twice the precession 
frequency and will be only a few percent (~ (3/2)z 2 ) of the overall amplitude. For 7 = 7r/4, 
Aq « 4i(/j,/D)(m/r) which is roughly 20% (initially) of the unmodulated amplitude in the 
previous case and is unmodulated through 0{i 2 ). The other three cases are modulated 
by a term with a frequency of twice the precession frequency, and another term with the 
precession frequency. For 7 = 7r/8, the dominant modulation is of 0(i) and varies with the 
precession frequency. In the other two cases (7 = i/2 and 7 = 2), both modulations are of 
0(i 2 ) and contain both a and 2a terms, which leads to their forms in Fig. |l[ Finally, let 
7 = i + 7r/4. Then Aq « 8i(fi/ D)(m/r)\ cos(a/2 — vr/4)| which corresponds with Fig. 6 of 
Ref. pi. 



V. RESULTS FOR SPECIFIC SYSTEMS 

In this section we will describe several specific cases involving a variety of masses and 
spins. For some cases, we have solved the equations of motion and precession numerically 
and used them to calculate the emitted gravitational waveform. In all the examples we will 
assume the binary's orbit has been circularized prior to entering the frequency bandwidth 
of a LIGO-type detector. 



A. Nonspinning Bodies 

The case in which neither body is spinning has been studied in previous papers 
Here we only wish to mention that if neither body is spinning, the orbital plane remains fixed 
in space, and the waveform's amplitude increases monotonically (apart from post-Newtonian 
modulations on orbital timescales, see Fig. |4j) as the binary inspirals. 
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B. Spins Perpendicular to the Orbital Plane 



If the spins of the bodies are aligned with the orbital angular momentum, the system 
evolves in a manner qualitatively similar to the case of nonspinning bodies. Since the spins 
and orbital angular momentum are aligned, none of them precesses so that the orbital plane 
remains fixed in space. The only effects the spins will have are a contribution to the orbital 
phase, and a correction to the amplitude of the waveform. 

The contribution to the orbital phase will be important as it affects the accumulated 
phase of the waveform. Since matched templates will be used to obtain information about 
the binary from the waveform, any effect which causes the phase to change by one cycle over 
the thousands in the bandwidth of the detector will be important. If the spin contribution 
to the phase can be separated from the other contributions, it could be used to make an 
accurate determination of the spins. However, preliminary studies by Cutler and Flanagan 



Tq suggest that the spin contribution cannot be separated cleanly from other contributions, 
but instead are strongly correlated with them, thus making the determination of the indi- 
vidual masses and spins more difficult. Table | compares the spin contribution to the orbital 
phase with other contributions for various binary systems. 

The leading-order correction to the waveform due to the spins is a full post-Newtonian 
order higher than the quadrupole part of the waveform. Thus it will be small until the late 
stages of the binary's inspiral. In Fig. |] we compare the spin contributions to the waveform 
with other contributions. 

C. One Spinning Body 

In the cases in which only one of the bodies is spinning, our numerical results agree with 



our analytic description of Sec. [TV|. As the binary inspirals, L and S precess about J, with 
both L and J shrinking, and L tipping away from J while S tips toward J (see Fig. |9|). The 
angle between L and S remains constant. J remains relatively fixed in direction unless L and 
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S are nearly antialigned and equal in magnitude, in which case the "transitional precession" 
described by Apostolatos et al |L9| occurs. 



1. Dependency of modulation on detector orientation and location 

The precession of the orbital plane causes the waveform amplitude to be modulated 
since the orientation between the orbital plane and the detector is changing. The form of 
this modulation depends on the orientation of the detector and its location with respect 
to the source. In Fig. [I], which was generated from numerical solutions, we show how 
the modulation changes for different orientations of a detector at a fixed location. Notice 
that the size and shape of the modulations vary greatly. These modulations are discussed 
qualitatively in Sec. |V[ 

2. Effects of higher order parts of the waveform 

So far, we have just examined how the precession of the orbital plane modulates the 
dominant quadrupole part of the waveform. The precession will also modulate the amplitude 
of the post-Newtonian corrections to the waveform. This is illustrated in Fig. |H| Recall 
that the different contributions to the waveform have different dependences on the orbital 
phase, so that the overall waveform amplitude is not the sum of the individual higher-order 
amplitudes. (Note that this is true even without spins.) 

D. Two Spinning Bodies 

We are not able to solve the general two-spinning-body problem analytically. We there- 
fore present numerical solutions of the equations of motion and precession. Most cases 
involving two spinning bodies are qualitatively similar to the case of one spinning body. 
The main difference lies in the fact that the total spin S is not constant as the spins precess, 
but rather oscillates between some S m i n and S max . This causes the angle i between L and 
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J to oscillate as L precesses about J in addition to its overall increase due to gravitational 
radiation damping. In most cases where spin effects will be important, these oscillations will 
be small. In each of the following examples, we will examine a pair of inspirals which have 
the same initial orientation of L and S, but in one case only one of the bodies is spinning, 
while in the other case both bodies are spinning. 

If one of the spins is much smaller than the other, it can be viewed as a perturbation of 
the simple precession involving the larger spin and the orbital angular momentum. Recall 
that the spin ratio will go as | Si | / 1 S2 1 = (xi/X2)("ii/'^2) 2 , so that if mi << m 2 and 
X2 is not small, then the spin of the smaller body will be much smaller than that of the 
larger body, and the above argument will hold. The overall shape of the modulations of the 
waveform amplitude is the same in both the one spin case and the two spin case. The main 
difference is that the precession frequency is slightly different in the two cases, which could 



lead to a noticeable effect on the phase of the waveform. Fig. [14] shows the difference in 
precession rates. 

If the masses of the two bodies are equal, Eq. (|2.5| ) implies that the only change in the 
evolution of S (and thus in |S|) compared to the one-spin case will be from the spin-spin 
coupling. This coupling is weak so that it can be viewed as a perturbation on the simple 
precession case. Fig. [15| illustrates this case. Notice that while the individual spins precess 
wildly and in effect exchange places, the total spin and orbital angular momentum remain 
approximately fixed relative to each other (re pa const.). This leads to the overall shape of 
the modulations of the waveform amplitude to be very similar for the one-spin case and the 
two-spin case. As in the previous example there is a difference in the precession frequency 
between the two cases, which could lead to a noticeable effect on the phase of the waveform. 

In some cases, however, the oscillations in i can be significant, and cannot be viewed as 
a perturbation of a simple precession case. Fig. [16] illustrates this case, for equal spins but 
different masses. Even though the spins are equal in magnitude, the different masses will 
cause them to precess at different rates (see Eq. ( |2.5| )). At some times the two spins will 
align and i will be at its largest value. At other times the two spins will almost cancel each 
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other, and L will be aligned with J. Thus the orbital plane will tilt back and forth as it 
precesses about J. This leads to substantial differences between the two spin case and the 



one spin case, as illustrated by the complex modulations in Fig. 16. 



Another case in which the second spin will be important is that of "transitional preces- 



sion" described by Apostolatos et al \ W\ . Recall that in this case L and S are almost cancel 
each other, so that any perturbations in S will cause noticable effects. See Ref. ]19|] for more 
details. 

VI. CONCLUSIONS 

In this paper we have examined the importance of spin-orbit and spin-spin effects on the 
inspiral of a coalescing binary system of compact objects and on the gravitational radiation 
emitted from such a system. The inclusion of spin effects makes the study of coalescing 
binary systems much more complicated because of the extra degrees of freedom for the 
orbit, and the extra parameters upon which the waveform will depend. On the other hand, 
if the effects due to the spins can be separated from other post-Newtonian effects, more 
information about the binary system can be extracted from the observed waveforms. 

The spins of the bodies have two major effects on the inspiral of the binary. As long 
as the spins are not perpendicular to the orbital plane, the orbital plane will precess, thus 
changing its orientation in space. In most cases, this precession will be a relatively simple 
precession of the orbital angular momentum L about the total angular momentum J. In 
some cases, however, the precession can be quite complicated (See Sec. In addition to 
the precession effect, the spins will contribute to the evolution of the orbital phase, in the 
same manner that other post-Newtonian terms contribute (See Sec. [IV] and Table |). The 



spin-orbit contribution can be of the same magnitude as the post-Newtonian and leading- 
order tail contributions to the orbital phase. The spin-spin contribution, on the other hand, 
is quite small. These contributions to the orbital phase will change the rate of inspiral. 
The spins will change the amplitude of the observed waveform in several ways. The major 
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effect the spins have is due to the precession of the orbital plane. This causes the orientation 
of the orbit to be changed with respect to the detector, which results in modulation of the 
amplitude of the waveform (See Sec. fVj, and Fig. [I]). The size and shape of the modulations 
are sensitive to the location of the detector with respect to the binary system and to the 
orientation of the detector arms. The spins also directly contribute to the amplitude of the 
waveform in the same manner that other higher post-Newtonian terms contribute (See Sec. 



IIJ and Fig. ||). The direct spin contribution to the amplitude may be difficult to detect 
because it is much smaller than the quadrupole term and the leading-order post-Newtonian 
terms until very late in the inspiral. The modulations due to precession, however, will be 
quite noticeable in many cases involving spinning bodies. 

The spins will also affect the phase of the observed waveform in several ways. Since the 
spins affect the evolution of the orbital phase, this will in turn affect the evolution of the 
phase of the gravitational waveform since it is related to the orbital phase. The phase will 
also be affected by the precession of the orbital plane since the point from which the orbital 
phase is measured is itself moving (See Sec. |TV| ). These two effects are independent of one 
another as one may be present without the other. These effects are very important, however, 
as the sensitivity with which the phase of the waveform can be measured is currently thought 
to be the best way of extracting information about the binary system. 

In general, the effects of the spins on the waveform amplitude will be small for the case 
of two coalescing neutron stars. One reason is that for the majority of the time in which 
the frequency of the gravitational waves from the binary are in the bandwidth of a LIGO- 
type detector, the orbital angular momentum will be much larger than the spin angular 
momentum. For example, over 95% of the gravitational wave cycles will occur between 
r = 174m and r = 37m for two 1.4M neutron stars, and as Fig. [5] shows, L is at least 5 
times larger than S. This means that the inclination angle i will be small throughout most, if 
not all, of the observed inspiral. This implies that the small inclination angle approximation 
will hold, so that the amplitude modulations will be on the order of i. For small enough i, 
one might be able to ignore the modulations to first order and treat the orbital plane to be 
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fixed perpendicular to J (which is relatively fixed) instead of L (which is precessing). Then 
the modulations can be treated as perturbations of 0(i). 

While the modulations of the amplitude of the waveform may be small, the spin effects 
on the phase of the waveform are significant for a coalescing binary system of neutron stars. 
While the spin-orbit contribution to the orbital phase is less than one percent of the total 
(see Table |), a neutron star with x — Sns/^ns as sman as 0.01 can cause the accumulated 
gravitational wave phase to change by a cycle from an equivalent system with no spins over 
the bandwidth of 10 Hz to 1000 Hz. A further contribution to the waveform's phase will come 
from the precession of the orbital plane. This contribution should vary at the precession 
frequency. Since L >> S for a binary neutron star coalescence, we may substitute L for J 
in Eq. (|2.14j) and integrate over the inspiral, finding that the binary precesses roughly 60-70 
times over the detector's bandwidth of 10 Hz to 1000 Hz. 

Spin effects can be very important for coalescing binary systems with a very massive black 
hole and a neutron star. Since the spins can dominate the orbital angular momentum, large 
precession angles are possible, leading to very large amplitude modulations. Furthermore, 
the spin-orbit contribution to the orbital phase is a larger percentage of the total. 

We have examined how the spins of the body affect the inspiral of the binary system 
and the gravitational radiation emitted therefrom. There are many interesting questions 
that remain to be examined. One is the inverse problem: given an observed waveform, how 
much information about the spins of the bodies can be extracted. A related question is how 
much will spin effects complicate the extraction of the masses and other information from 



the waveform. Preliminary studies by Cutler and Flanagan []17| indicate that the spin-orbit 
contribution to the orbital phase will be difficult to separate from the post-Newtonian con- 
tribution which is used to determine the mass ratio because of their similar dependences 
on (m/r) or (moo) in the evolution of the phase (see Eq. |4.16| ). Further studies are needed 
to determine whether Newtonian templates without explicit spin contributions to the wave- 
forms can be used on systems which have spins without significant loss of signal-to-noise 
ratio. If not, then the number of templates which will be needed to study coalescing binary 
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systems of compact objects will increase significantly. 
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When examining systems containing spinning bodies, it is important to note that the 
form of the spin-orbit acceleration &so is not unique, but rather depends on a "spin supple- 
mentary condition" (SSC). For example, three different forms of &so 



APPENDIX A: SPIN SUPPLEMENTARY CONDITIONS 
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are given by the three different SSC's 



S^uau = 0, (A2a) 
2S Ai0 + S Ai y A = 0, (A2b) 



S% = 0, (A2c) 

respectively, where u A is the four-velocity of the center-of-mass world line X A of body A, 
and 

S A V = 2 / (a> - X A ^)r u]0 d 3 x, (A3) 



.4 

where t^ u denotes the stress-energy tensor of matter plus gravitational fields satisfying 
rfJ,u ,u — 0, and square brackets around indices denote antisymmetrization. Note that the 
spin vector S of each body is defined by S A = ^ijkS A - 

Barker and O'Connell ]27| showed that these different forms of slso are equivalent if 
one takes into account that the different SSC's are related to different locations Xa of the 
center of mass of each body. Furthermore they found transformations from the center-of- 
mass definition given by one SSC to that of another. For example 

X\ [n) — > Xf ] l^-(v A x S A y , (A4) 

This shift in the center-of-mass world line is of post-Newtonian order, so it can be neglected 
at lowest order. 

In Sec. [TT], we chose to use the form of &so given by the first SSC (|A2a[J , since it is 



covariant. In doing so, we have chosen a center-of-mass definition through the SSC, and we 
must insure that any future center-of-mass definitions we use are consistent with this choice. 



In Sec. |T|, we use an integral definition of the center of mass of each body, 

d 3 x, (A5) 



X ' A = I %lp * (x) 

A 
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to evaluate the BDI multipoles, where 



m A = l p (xj 

A 



d 3 x, 



(A6) 



where v\ = v 1 — v\, v\ = dx A /dt, and U A is the Newtonian potential produced by the A-th 
body itself. This is a "natural" post-Newtonian definition since it uses the total mass, kinetic 
energy, internal energy, and gravitational potential energy of each body as the weighting 
factor. However, it turns out that this definition of the center of mass is not related to the 
first SSC which we chose for our equations of motion, but rather to the second SSC (|A2b| ) . 
This can be seen by the following argument. 

It is straightforward to show that the three SSC's can be rewritten in the form 



£>A ~ ^^a v a ~ 0' 



(A7) 



where k = 1 for the first SSC, k = 1/2 for the second, and k = for the third. From Eq. 
we see that 



J A 



(A8) 



since the integration is done at constant time. It is straightforward to show that, to post- 
Newtonian order, 



7 



r uu = (1 + 417) T uu -— |Wr, 

8n 



(A9) 



Substituting Eq. (|A9|) into Eq. (|A8|) and integrating the |VC/| terms by parts, we obtain 



1 + i v 2 _ l _u + n 
2 2 



d 3 x. 



(AlO) 



Using our integral definition of the mass and center of mass of body A (|A5|) , we obtain 
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Thus our integral definition of the center of mass corresponds to the SSC definition (x A 



\v A S 1 ^ i.e. the second SSC. Imposing the general SSC 



on Eq. flSIID , 



we see that, to post-Newtonian order, 



X A + 



2Jfc- 1 
2mA 



Q«3 



A v Ai 



(A12) 



Thus in order that the center of mass used in the STF moments be consistent with that used 
in the equations of motion, we must transform from x\ = X A ir * to X A in the moments by 
using Eq. (|A|) (See Sec. g). 

This is also consistent with Brumberg [IT who derived the equations of motion from our 
integral definition of the center of mass, and found the form of slso to be that of the second 
SSC without ever having mentioned the concept of an SSC. 

Another check of the consistency of our argument is seen by evaluating the BDI mass 



dipole moment [23,24], given by 



E 



1 9 

1 + -v 2 



lu + u 



d 6 x 



(A13) 



Evaluating the integral in the manner used to evaluate the mass quadrupole (in Sec. [TTJ), 
we obtain 



A V 



1 2 1 m B 

1 1 B^A T AB 



1 



+ - (v A x S A ) 



(A14) 



where we have dropped the post-Newtonian terms as they do not affect the calculation. We 
then use Eq. (|A4|) to shift the center of mass x A from that defined by the second SSC X l A ^ 



to that of our equations of motion, X A . The result is 



Ji = Ei m A x A 

A I 



1 9 

i + -vi 



1 x - % 

2 £ A T AB 



+ (v A x S A ) 



(A15) 



where x l A now refers to X\ . Blanchet and Damour [23] have proved on general grounds 
that 
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dt 2 ' 



(A16) 



Taking two time derivatives of Eq. ( |A15| ) and substituting our equations of motion where 
appropriate, we verify that Eq. ( |A16| ) does hold. (Alternatively, we could have taken two 
time derivatives of the mass dipole before shifting the center of mass of each body and seen 
that SigQ is needed to insure that Eq. ( |A16| ) holds.) 



APPENDIX B: POST-NEWTONIAN CORRECTIONS TO THE WAVEFORM 

FOR A CIRCULAR ORBIT 



If we define $ as the orbital phase with respect to the line of ascending nodes (the point 
at which the orbit crosses the x-y plane from below), then for a circular orbit the polarization 
states will be given by 



h+ M\[ Q++ ^ Q+ {^\ PQ jvi 



where, for simplicity we include only the post-Newtonian terms through PQ + , and 



(Bl) 



Q- 



-2 [C+ cos 2$ + S + sin 2$] 



(B2a) 



1 5m 



P U5 Q + = — [9 (aS + + bC + ) cos 3$ + 9 (bS + - aC + ) sin 3$ 
4 m 

+ (3aS + - 3bC + - 2bK + ) cos $ - (3bS + + 3aC + - 2aK + ) sin$], 



(B2b) 



PQ + = _(l - 3?7) {[(a 2 - b 2 )C + - 2abS + ] cos 4$ + [(a 2 - b 2 )S + + 2abC+] sin 4$} 
+L>CV cos $ + sin $ + -([4(1 -3r^)(a 2 + 6 2 ) + (19-3^)1 Q+ 

6 L L J 

-4(1 -377) [(a 2 -6 2 )cos2$ + 2a6sin2$] K + }, (B2c) 
and similarly for h x where + is replaced by x, and where 

ri 1 2 r\ ( ■ 2 2- 2\,1/ 2 • • 2 

G + = 2 cos B I sin a — cos «cos a J + | (cos ism a 



(B3a) 



cos 2 a) — \ sin 2 B sin 2 i — | sin 2B sin 2i cos a, 
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S + = | ( 1 + cos 2 0J cos i sin 2a + | sin 20 sin i sin a, (B3b) 

TS 1 2 d ( ■ 2 i 2- 2 \ 1/ 2 • • 2 

A + = icos B Ism a + cos 2 cos a J — ^(cos 2 sin a 
+ cos 2 a + | sin 2 sin 2 z + ~ sin 20 sin 2i cos a, 

(B3c) 



DC+ = [A y sin a cos + d cos a] , (B3d) 

rrr 

DS + = \cA y — d cos i sin a] , (B3e) 

C x = — | cos sin 2a (l + cos 2 ij — | sin sin 2z sin a, (B3f ) 

= — cos cos i cos 2a — sin sin i cos a, (B3g) 

fT x = — ~ cos sin 2a sin 2 z + ~ sin sin 2z sin a, (B3h) 

L>C X = — - \A y cos a; — d cos sin a] , (B3i) 

D^x = -A? [-A^coszsina + cd] , (B3j) 
m 2 

a = — sin sin a, (B3k) 

6 = cos sini — sin cosi cos a, (B31) 

c = cos cosi cos a + sini sin 0, (B3m) 

d = A 2 sin© - A x cos0, (B3n) 



where (i, a) are the spherical coordinates describing the direction of L N (see Fig. §). 
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From the above equations, it can be seen that the signal in the detector can be written 
as a series 

2(N+1) 

h= i C n cos n$ + S n sin n$} , (B4) 

71=1 



which can be rewritten as 



2(JV+1) 

h= A n cos[n$ - 5 n ], (B5) 

n=l 

where 

A n = {c 2 n + Sl} l/ \ (B6a) 

5 n = tan- 1 (S n /C n ) , (B6b) 

and N is the final post-Newtonian order beyond the quadrupole term to which the waveform 
is calculated. For example, if we look at only the quadrupole contribution to the waveform, 
then only the n = 2 terms will be nonzero and 

C 2 = Jt (™) [ C+ F + + C X F X ] , (B7a) 



D \ r 



S 2 = (™) [S + F + + S X F X ] . (B7b) 



From Eqs. (|B2j ) we see that the post-Newtonian corrections to the waveform depend 
on different harmonics of the orbital frequency than the quadrupole term. Eq. ( |B5|) splits 
the signal into the terms depending on the different harmonics. Each term will have its 
amplitude modulated by the precession of the orbital plane. 

APPENDIX C: A COMPARISON WITH TEST MASS CALCULATIONS 

Recently, Poisson |34| has calculated the gravitational waveform and energy loss due 



to gravitational radiation from a particle in a circular, equatorial orbit around a slowly 
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rotating black hole. Poisson's calculations were done using the Teukolsky perturbation 
formalism, which is completely different from the post-Newtonian calculations we used. In 
this appendix, we show that our results for the waveform and energy loss agree with those 
of Poisson in the appropriate limit. Since we are only interested in the spin-orbit terms, we 
will neglect all the other post-Newtonian terms. 
Poisson gives the energy loss as 



d_E = _Z2f }L] (// ^„ 1 , ! 
at 5 V m . 



1 - -^XBH{muJ) 



(CI) 



and the spin-orbit contributions to the waveform as 



h + = 2 -§(mu)^\Q + - XBH 



sin sin $(mu) 2 ^ 3 



- - ( 1 + cos 2 ) cos 2$(mco>) 
3 



(C2a) 



2 At, 
~~D 



[moo 



,2/3 



Qx + Xbh 



sin cos cos $(mw) 2 ^ 3 



- cos sin2$(mc(/ 
3 V ' 



where we have rewritten the equations using our notation. 

Combining Eqs. ( [4.10|) and ( |4.13| ) we find the energy loss to be 



(C2b) 



dE 



32 



m 



4, 



Xi(L N ■ §i) 



i=l,2 



x 11—^ + 5?7 (mw) 



(C3) 



In the test-mass limit that we are considering, ni\ = fi, rri2 = m (with \i << m), 
rj = fi/m, xi — 0, X2 = Xbh, Sm = —m, s 2 = Ln, and A = xs^m 2 LN- Substituting 
these values into Eq. ( |C3D we see that our expression for the energy loss agrees with that of 
Poisson in the appropriate limit. 

Combining Eqs. (O) and (|4.13 ), and taking the above test mass limit of the spin 
contributions yields 

\2/3 



+ ^XBH{muj) 



n l n 3 - X\ d 



(C4) 
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Substituting n = cos $x + sin $y and A = — sin $x + cos $y, into the above expression and 
inserting the results into Eqs. (|4.23|) we obtain a result which matches that of Poisson. 
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FIGURES 

FIG. 1. Amplitude modulation of gravitational waveforms by spin-induced orbital precession, 
plotted against time to coalescence. System consists of a nonspinning 1M & neutron star and a 
maximally spinning black hole of 1QM®. Spin and orbital angular momentum vectors are initially 
misaligned by 11.3°. Initial orbital inclination relative to J is i. The angle 7 represents the 
orientation of the detector relative to J, with the detector located on the x-axis (see Fig. ^) such 
that the source is directly overhead. Curves show envelope of the quadrupole waveform for various 
detector orientations (The curves for the cases 7 = 0, i/2, and i would lie on top of one another, 
so the first two have been shifted upward for ease of presentation). Gravitational- wave frequency 
runs from 10 Hz on the right to 300 Hz on the left. 

FIG. 2. The source coordinate system. The total angular momentum J initially lies along the 
z-axis. The detector is located in the x-z plane. The spherical angles (i, a) define the direction of 
the Newtonian angular momentum Ln which is perpendicular to the orbital plane. In terms of 
celestial mechanics, the angle of ascending nodes is a + ir/2. 

FIG. 3. A comparison of the magnitudes of the orbital angular momentum and spin angular 
momentum as the binary inspirals. (a) The equal mass case (assuming the bodies are maximally 
spinning), (b) The case of a 10:1 mass ratio. 

FIG. 4. Gravitational wavefrom plotted against orbital phase for a 10:1.4 mass-ratio system. 
The smaller body's spin is aligned with the orbital angular momentum L, while the larger body's 
spin is tilted by an angle of 30° with respect to L. Plotted is (D/2(i)h+ for an observer at = 90°. 
Plots begin at an orbital separation of 15m and terminate at 10m. (a) The total waveform, 
(b) The quadrupole contribution to the waveform, (c) The first higher-order post-Newtonian 
correction ((^(e 1 / 2 ) beyond the quadrupole). (d) The next post-Newtonian correction, (e) The 
leading-order spin-orbit contribution to the waveform, (f ) The leading-order spin-spin contribution 
to the waveform (note the different scale) . Notice the modulation due to the precession of the orbital 
plane. 
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FIG. 5. The wobble of the orbital plane during simple precession (in the absence of gravi- 
tational-radiation damping). Plotted are the x- and y-components of the unit vector Ln which 
defines the orbital plane. The total angular momentum J is directed out of the page. If the orbital 
plane were not wobbling as it precessed, the plot would be a circle whose radius depends on the 
inclination of the orbital plane with respect to J. 

FIG. 6. The complicated precession of the orbital angular momentum L (in the absence of 
gravitational-radiation damping) for a 2:1 mass-ratio system with equal spins which are initially 
aligned in the orbital plane. Although the spins are equal, they precess at different rates so that 
at a later time they almost cancel one another. This leads to the tilting of the orbital plane as it 



precesses about J. See Fig. 10 for the effects of gravitational radiation on the precession. 



FIG. 7. The energy lost due to gravitational radiation plotted against orbital phase for a 10:1.4 
mass-ratio system (with spins initially aligned as in Fig. ^]). Plots begin at an orbital separation 
of 15m and terminate at 10m. Plotted are the lowest-order Newtonian (N), post-Newtonian (PN), 
spin-orbit (SO), and spin-spin (SS) contributions to the energy loss. 

FIG. 8. The linear momentum ejected due to gravitational radiation plotted against orbital 
phase for a 10:1.4 mass-ratio system (with spins initially aligned as in Fig. ||). Plots begin at an 
orbital separation of 15m and terminate at 10m. 

FIG. 9. The simple precession of the orbital angular momentum L in the presence of gravita- 
tional radiation damping for an inspiraling system with a 2:1 mass-ratio with only the larger body 
spinning (The spin is perpendicular to L). Notice that the inclination of the orbital plane with 
respect to J Q increases as the binary inspirals (i = sin _1 [(L 2 + L 2 ) 1 / 2 ]). 

FIG. 10. The complicated precession of the orbital angular momentum L in the presence of 
gravitational radiation damping for an inspiraling system with a 2:1 mass-ratio and equal spins 
which are initially aligned, and perpendicular to L. Notice that the maximum inclination of the 
orbital plane increases as the binary inspirals. 
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FIG. 11. The detector coordinate system. The detectors arms lie along the x- and y-axes. The 
angles (6, <p) define the location of the source, and the polarization angle ip describes the orientation 
of the polarization axes with respect to the projection of the detectors arms on the sky. 

FIG. 12. The region of validity for the small inclination approximation for an equal mass 
system. The figure plots the minimum separation r/m for which the precession angle i is guaranteed 
to be less than i ma x as a function of the angle between L and S for several values of the average 
spin parameter x °f the two bodies. For most neutron star models x ^ 0.7. Note that at large 
angles, "transitional precession" may occur. 

FIG. 13. The amplitude modulations of the higher-order terms of the gravitational wave- 
form. A\ and ^3 are of (post) 1//2 -Newtonian order beyond the quadrupole term, while is 
(post) 1 -Newtonian order. The quadrupole envelope increase from 0.06 to 0.36 during the portion 
of the inspiral which is plotted. Note that since the different terms have different dependences on 
the orbital phase, the overall waveform envelope will not be the sum of the individual envelopes. 
The binary system is the same as in Fig. [l], but for an observer at = 7r/4, with 7 = 0. 

FIG. 14. The rate of precession for the two-spin case and the corresponding one-spin case for a 
10:1 mass-ratio system with the total spin S at an angle of 7r/4 with the orbital angular momentum 
L. In the two-spin case, the smaller body's spin is aligned with L. In both cases the bodies are 
maximally spinning. 

FIG. 15. An equal mass case in which the spins of the bodies are equal in magnitude 
(xi = X2 = \/2/2), with one spin aligned with L and the other perpendicular to it. (a) The 
modulation of the quadrupole waveform envelope, (b) The components of the total angular mo- 
mentum J. Notice that J remains relatively fixed in direction as the binary inspirals due to 
gravitational radiation emission. Also plotted is the angle between L and S, which oscillates about 
a relatively constant value, (c) The components of Si. (d) The components of S2. Notice that 
while the individual spins undergo substantial precessions, the total spin precesses in a fairly simple 
manner. This causes the two-spin case to be similar to the corresponding one-spin case. 
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FIG. 16. The modulation of the quadrupole waveform envelope for a 2:1 mass-ratio system with 
the total spin S initially perpendicular to L. (a) One spinning body, (b) Two spinning bodies with 
spins equal such that the smaller body's spin is maximal (x = 1). Note the substantial differences 
in the modulations between the two cases. This is a result of the complicated precession when both 
bodies are spinning. See Figs. |9] and 10 for a depiction of the precession of the orbital angular 
momentum for these cases. 
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TABLES 

TABLE I. The contributions to the number of orbits for an inspir ailing binary system which 
is in a LIGO-type detector's bandwidth, from quadrupole (Quad), post-Newtonian (PN), tail, 
spin-orbit (SO), and spin-spin (SS) terms. Masses are in units of M Q . In the first table, the 
observed inspiral begins when the gravitational wave frequency enters the bandwidth at the seismic 
cutoff around 10 Hz, and is cut off when either the quadrupole radiation of the binary has left 
the detector's bandwidth (at around 1000 Hz), or when the binary reaches its innermost stable 
circular orbit at which point the binary will plunge to a single object in one or two orbits. In the 
second table, we list the contribution to the number of orbits in the narrower bandwidth of 40 Hz 
to 100 Hz. The spin contributions shown are the maximum contributions assuming the bodies are 
maximally spinning (% = Xmax)- The magnitude and sign of the spin contributions depend on the 
specific orientations of the spins with respect to the orbital plane. 
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(f) Spin-spin P 2 N xlOO 
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Angle Between L N and S (deg) 



